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.it  reflected  enerqy/incident  enerqy)  computed  from  spectra  ratios  vary  stronq 
with  incidence  angle.  At  wavelengths  equal  to  twice  the  step  height.,  the 
fraction  of  incident  energy  scattered  into  body  waves  ranges  from  more  than 
'Or  at  normal  incidence  to  about  zero  at  near-grazing  incidence.  At  each 
angle,  transmission  coefifcients  vary  strongly  with  frequency.  Because  of 
frequency-dependent  phase  shifts,  the  transmitted  and  reflected  waves  are 
d  istorted . 

The  effect  of  the  steps  on  the  propagation  of  Rayleigh  waves  is 
lemonst rated  by  convolving  synthetic  dispersed  wave  trains  with  the  impulse 
iespon.se  of  the  scale  models.  The  ocean-continent  margin  of  the  western 
United  States  is  modeled  as  a  60°  ramp  scaled  to  60  km  height.  The  Tibetan 
Plateau  is  modeled  as  a  broad  mesa  scaled  to  40  km  height.  In  both  models 
the  azimuthal  dependence  of  transmitted  Rayleigh  waves  is  similar  to  that 
observed  at  WWSSN  stations  for  Rayleigh  waves  crossing  the  modeled 
terrestrial  structures. 


ABSTRACT 


The  effects  of  topographic  features  on  Rayleigh  wave  propagation  and  scatter¬ 
ing  are  investigated  in  the  laboratory  using  three-dimensional  ultrasonic 
models.  Starting  from  simple  steps,  different  topographic  features  are  modeled. 
The  effects  of  these  features  on  Rayleigh  wave  transmission  and  scattering  arc 
examined  as  a  function  of  wavelength  nd  as  a  function  of  angle  of  incidence. 
In  general,  backscattered  or  reflected  Rayleigh  waves  are  small  compared  to 
transmitted  waves.  A  significant  fraction  of  the  Rayleigh  wave  energy  is  scat¬ 
tered  into  body  waves.  Transmission  and  reflection  coefficients  (transmitted  or 
reflected  energy/incident  energy)  computed  from  spectral  ratios  vary  strongly 
with  incidence  angle.  At  wavelengths  equal  to  twice  the  step  height,  the  frac¬ 
tion  of  incident  energy  scattered  into  body  waves  ranges  from  more  than  90%  at 
normal  incidence  to  about  zero  at  near-grazing  incidence.  At  each  angle, 
transmission  coefficients  vary  strongly  with  frequency.  Because  of  frequency- 
dependent  phase  shifts,  the  transmitted  and  reflected  waves  are  distorted. 

The  effect  of  the  steps  on  the  propagation  of  Rayleigh  wave?' is  demon¬ 
strated  by  convolving  synthetic  dispersed  wave  trains  with  the  impulse  response 
of  the  scale  models.  This  is  done  using  the  transmission  response  function  at 
the  appropriate  angle  and  from  the  simple  model  described  above  and  convolv¬ 
ing  with  the  input,  wave  function.  The  ocean-continent  margin  of  the  western 
United  States  is  modeled  as  a  60®  ramp  scaled  to  60  km  height.  The  Tibetan 
Plateau  is  modeled  as  a  broad  mesa  scaled  to  40  km  height.  In  both  models  the 
azimuthal  dependence  of  transmitted  Rayleigh  waves  is  similar  to  that  observed 
at  WWSSN  stations  for  Rayleigh  waves  crossing  the  modeled  terrestrial  struc¬ 
tures.  The  actual  physical  modelling  of  the  plateaus  is  now  underway  using  sim¬ 
ple  shapes  (such  as  circular  mesas)  at  first. 
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INTRODUCTION 


Rayleigh  wave  propagation  in  strongly-varying  structures  is  a  complex  pro¬ 
cess  well  suited  to  analysis  by  ultrasonic  modeling  methods.  Until  recently, 
ultrasonic  modeling  of  surface  wave  transmission  and  scattering  has  been  lim¬ 
ited  to  two  dimensions.  In  this  paper  we  report  on  a  series  of  ultrasonic  model¬ 
ing  experiments  investigating  the  propagation  and  scattering  of  Rayleigh  waves 
by  simple  three-dimensional  topographic  features. 

Rayleigh  waves  propagating  in  two  dimensions  through  simple  steps  or 
around  corners  and  in  wedges  have  been  studied  experimentally  by  a  number  of 
investigators  (dc  Bremacckcr,  1958;  Kato  and  Takagi,  1956:  Lewis  and  Dally, 
1970;  Knopoff  and  Gangi,  I960;  Pilant  et  al.,  1964;  Martell  et  at.,  1977;  and  Nath- 
man,  1980}.  Detailed  'snapshots'  of  the  complete  stress  field  of  Rayleigh  waves 
scattered  in  wedges  have  been  obtained  using  photoelastic  methods  (Lewis  and 
Dally,  1970). 

Development  of  theoretical  and  numerical  methods  of  synthesizing  Ray¬ 
leigh  waver  scattered  from  simple  structures  has  typically  lagged  behind  exper¬ 
imental  results.  Rayleigh  wave  transmission  and  reflection  coefficients  across 
vertical  boundaries  in  two  dimensions  have  been  calculated  using  approximate 
variational  methods  (e.g.  McGarr  and  Alsop,  1967).  Recently,  synthesis  of  com¬ 
plete  seismograms  of  Rayleigh  waves  reflected  from,  and  transmitted  through,  a 
vertical  step  has  been  performed  using  the  finite-difference  method  (Toksoz 
198.!;  Fuyuki  and  Nakano,  1984).  In  three  dimensions,  numerical  computations 
arc  limited 

Approximate  methods  to  estimate  transmission  and  reflection  coefficients 
of  Raylt.igl  waves  incident  at  steep  angles  on  a  vertical  boundary  have,  how¬ 
ever,  been  developed  (Maliehewsky,  1976;  Chen  and  Alsop,  1979). 


Scattering  of  Rayleigh  waves  in  three  dimensions  can  be  investigated  in  the 
laboratory  using  ultrasonic  methods.  Rayleigh  wave  propagation  in  three- 
dimensions  has  been  studied  using  realistic  three-dimensional  scale  models  of 
surface  topography  (e.g.  Toksoz,  1983).  In  such  models  even  a  simple  input 
pulse  is  severely  distorted  by  the  complex  structure.  Subtle  effects  of  topogra¬ 
phy  on  Rayleigh  wave  propagation  and  scattering  are  thus  easily  obscured  by 
the  complexity  of  the  model.  An  important  step  toward  an  understanding  of 
surface  wave  scattering  in  three-dimensions  is  the  investigation  of  propagation 
in  simple  models  of  canonical  form. 

In  this  paper,  we  investigate  Rayleigh  wave  scattering  from  step  and  ramp- 
type  topographic  discontinuities  and  study  their  implications  for  surface  wave 
propagation  across  major  features  such  as  ocean-continent  boundaries  and 


EXPERIMENTAL  METHOD 


••j 

The  model  used  in  this  study  consists  of  a  solid  block  (200x200x100  rnirr) 
of  aluminum.  Steps  of  various  height  and  geometry  are  milled  into  one  of  the 
large  fanes  (Figure  1). 

The  ultrasonic  source  and  receiver  each  consist  of  a  2.25  MHz  P  wave 
transducer  coupled  to  a  lueite  wedge.  The  geometry  of  the  wedges  is  such  that 
P  wave's  generated  by  the  transducer  excite  Rayleigh  waves  in  the  aluminum 
with  high  efficiency.  Similarly,  the  transducer  is  an  efficient  Rayleigh  wave  sen¬ 
sor.  The  radiation  pattern  of  the  transducer- wedge  combination  is  highly 
directional;  at.  an  angle  of  10°  from  the  axis  of  the  central  lobe,  peak  energy  is 
reduced  to  less  than  25%  the  axial  peak  energy.  The  combination  of  these 
features  makes  the  transducer- wedge  combination  especially  useful  when  work¬ 
ing  with  a  low-loss  medium  such  as  aluminum;  attenuation  can  be  neglected  in 
the  analysts  of  Rayleigh  wave  propagation,  and  the  observed  Rayleigh  waves  are 
ure-ont animated  by  body  waves  multiply  reflected  at  the  edges  of  the  model. 

Peak  energy  of  Rayleigh  waves  traveling  on  the  aluminum  block  ( =  (5.4 
km/ sec;  P,  =  3.2  km/sec;  =  3.0  km/sec)  lies  in  the  range  0.2— 2.0  MHz 
(A=  1.5—1  5  rnrn).  In  this  frequency  band  the  spectrum  is  repeatably  obtained 
from  a  given  model  and  source-receiver  configuration.  We  therefore  take  this 
as  the  usable  bandwidth  of  the  modeling  experiments.  Because  the  model  is 
essentially  a  half-space,  only  fundamental- mode  Rayleigh  waves  are  observed. 
Since  a  half-space  is  non-dispersive  to  Rayleigh  waves,  the  half-space  Rayleigh 
wave  r  epresents  the  impulse  response  of  t  he  transducers  and  recording  system. 

When  a  Rayleigh  wave  encounters  a  step  some  of  the  incident  energy  is 
ti  aiiMiiitled  through  the  step,  some  is  r  eflected  from  the  step,  and  some  is  ron- 
v"rtcd  to  body  waves  (de  Hremaecker,  1953;  KnopofT  and  Gangi,  19(50;  Martel  at 
rd  1 977,  Ink  so/.,  19H3).  In  general  the  effc  leney  of  transmission  through  m 
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reflection  from  a  step  depends  on  the  wavelength  relative  to  the  step  height 
(e.g.  McGarr  and  Alsop,  1967).  In  order  to  facilitate  comparison  of  models  with 
different  step  heights,  therefore,  wavelengths  are  normalized  to  the  step  height. 
The  "normalized  frequency",  / ,  is  obtained  from  the  scaling  relation 


1  c  X 

where  f  is  the  frequency,  e  is  the  phase  velocity,  X  is  the  wavelength,  and  h  is 
the  step  height.  For  a  given  step  height.,  the  range  of  usable  normalized  fre¬ 
quencies  is  limited  by  the  usable  bandwidth  of  the  modeling  experiments.  This 
is  ari  important,  consideration  when  the  model  seismograms  are  scaled  to 
correspond  to  real  Flartti  seismograms.  Transmission  and  reflection  coefficients 
calculated  from  Rayleigh  waves  traveling  across  steps  of  various  heights  indi¬ 
cate  that  the  model  scaling  relation  (Equation  l)  is  generally  appropriate  over  a 
wide  range  of  frequencies  and  step  heights. 


Transmitted  Rayleigh  waves  are  recorded  by  keeping  source  and  receiver 
at  opposite  points  on  the  diameter  of  a  circle  centered  on  the  step  (see  Figure 
I)  Reflected  waves  are  observed  by  placing  source  and  receiver  at  equal  angles 
from  t  tie  normal  to  the  strike  of  the  step.  Incidence  angle  is  measured  from  the 
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MODELING  RESULTS 


VERTICAL  STEP  MODEL 

Figure  3  shows  seismograms  of  Rayleigh  waves  transmitt ed  arrow  a  hmrn 
vertical  step  at  incidence  angles  ranging  from  15°  through  90°  (normal 
incidence).  The  small  signal  arriving  about  1  /isee  before  the  Rayleigh  wave  is 
an  S  wave  excited  at  the  step.  At  a  given  incidence  angle  the  waveforms  for  111' 
case1  of  up-  and  down-step  transmission  are  essentially  identical.  Not.tr  that  not 
only  peak  amplitudes  decrease  by  a  factor  of  two  to  four  or  more,  but  also  the 
waveform  changes  strongly  as  a  function  of  angle  of  incidence.  At  normal 
incidence  (0=90°)  transmission  coefficients  (Figure  3)  agree  well  with  those 
obtained  for-  a  two-dimensional  step  using  theoretical  (numerical)  calculations 
(Drake.  1973;  Martel  al  ,  1977;  Fuyuki  and  Nakano,  1984). 

At  normal  incidence  t.he  fraction  of  transmitted  energy  (EV)  seldom 
exceeds  30%  within  the  usable  frequency  band.  This  represents  a  loss  of  over 
H0%  of  the  incident  Rayleigh  wave  energy  to  reflected  surface  waves  and  con 
verted  body  waves. 

Reflected  Rayleigh  wave  seismograms  are  complicated  by  the  presence  of 
double  reflections,  separated  by  about  3  jj. see  (Figure  4).  In  Figure  4A,  the  first 
reflection  corresponds  to  the  Rayleigh  wave  reflected  from  the  bottom  edge  of 
the  step,  the  second  to  t.he  Rayleigh  wave  reflected  from  the  upper  edge  after 
propagating  up  t.hc  face  of  the  step.  These  are  further  complicated  by  multiple 
scattering  and  Rayleigh  to  shear  scattering.  Reflections  from  a  down-step  (Fig- 
ur>  1H)  arc  much  larger  than  those  from  an  up-step  (Figure  4A).  Also,  shapes  of 
dnw->-«tep  reflected  waveforms  are  different  from  those  of  up-st.ep  They  do  not 
exhibit  t.he  strong  incidence  angle  dependence.  Similar  differences  in  up- step, 
down  •  vt  cp  s—att  ering  of  Rayleigh  waves  have  been  observed  in  two-dimensional 


modeling  experiments  (e.g.  Nathman,  1990). 

In  general  the  reflected  energy  is  very  small,  particularly  at  normal  or 
near-normal  incidence  (Figure  (S).  At  most  angles  and  frequencies,  reflected 
energy  accounts  for  no  more  than  10-20“-?  of  the  incident  energy  As  the 
incidence  angle  becomes  more  grazing  the  reflected  energy  increases  rapidly, 
especially  for  wavelengths  approximately  equal  to  three  times  the  step  height. 
if  'i> .'FI 

It  ;s  clear  from  Figures  3  and  f>  that  the  efficiency  of  body  wave  scattering 
.  d  rongly  dependent  on  the  angle  of  incidence.  In  the  3  mm  step  model,  at 
wavcii  ngths  equal  to  one-half  the  step  height  (the  first  maximum  in  Figure  3) 
•.O'.  ;*  >;■ "  of  t  tie  incident,  energy  is  scattered  to  body  waves.  As  the  incidence 
. .  to  near  ()r  (gr  azing  incidence),  this  figure  drops  to  about  zero. 

v '  1  i !■;/ 

i'o  invent  ig.it  e  Itie  effect,  of  a  more  gradual  change  of  elevation,  the  3mm 
a  i  •:  tap<  red  to  a  slope  of  80°  from  the  horizontal  (see  Figure  1).  Distortion 
>‘  no’  n  transmitted  and  reflected  waveforms  is  strong,  particularly  at  normal 
and  "",:r  normal  incidence:  (Figure  6). 

Fncrgy  transmission  eot  ffleient  s  for  this  model  (Figure  7a)  are  generally 
■r  id  an  t  hose  for  the  step  mode)  (Figure  3b)  at.  frequencies  below  /  ^  0.f>. 

.  t  >•  "rgy  minimum  near  f  *0.9  -it  normal  incidence  is  considerably  smaller 
Ida  i  >  •••orr»,>p»nr|ing  minimum  fo*  the  step  model.  Reflected  energy  (Figure 
;t  d  .o.-ader  ably  smaller  than  for  I  tie  ease  of  (he  vertical  step  (Figure  4a). 

.  f  i  do-coir  reflect.  >n  :s.  howi-ver,  -if  parent  at  all  angles.  Interferenc  between 
ft  «•  two  reflections  and  t  he  scattered  body  waves  contribute  to  the  •  •omplcxd  y 
-  >1  t  h>-  seismograms 
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Figure  ID.  The  region  in  h-A  space  in  which  the  ultrasonic  model  results  from 
the  drum  models  may  be  reliably  used.  Units  of  h  and  A  are  arbitrary.  In  the 
shaded  areas  repeatability  and  consistency  of  a  given  experiment  is  not 
assiip  d.  Th  c  usable  and  unusable  regions  are  demarcated  by  the  maximum  and 
mimmum  usable  normalized  frequencies.  The  box  shows  the  region  allowed  for 
modeling  Rayleigh  waves  in  the  period  range  lb-40  see  (hO-lhO  km  wavelengt  h). 
Thus  step  heights  in  the  range  30-100  km  may  be  modeled  successfully  over 
tins  period  range. 

Figur  e  11.  A.  WWSSN  vertical  component  recordings  of  the  Rayleigh  waves  from 
an  earthquake  in  the  Loyalty  Islands  (9/6/91.  1  lh03m40.9S,  31.4H9r  IS,  109.000 
L.  <1  ept  h  —  3 1  km,  rn,  =f>  9,  Ms=6.3).  The  angle  at  which  the  great-circle  crosses 
the  continental  margin  is  shown  next  to  the  continental  stations.  B.  Model  syn¬ 
thetic  seismograms  using  the  3mm,  60°  ramp  scaled  to  60  krn  height. 


Figure  13.  Rayleigh  wave  paths  across  the  Tibetan  Plateau  from  an  explosion  .it 
the  Lop  Nor  test  site  (14  October,  1970).  Reproduced  from  Bird  (1976).  Inset: 
lnr-ce-dimensional  model  of  the  Tibetan  Plateau  used  to  generate  model- 
s> nt Let ic  seismograms  for  the  paths  Lop  Nor-NDI  and  Lop  Nor-SHL. 


Ugure  H  Comparison  of  real  seismograms  (solid  lines)  with  model-synthetic 
s*  . -mograms  (dashed  lines)  using  the  lens  model  of  Figure  13.  All  seismograms 
base  been  low-pass  filtered  at  10  second  period. 


FIGURE  CAPTIONS 


Figure  1.  Geometry  of  the  models  used  in  this  study.  Source  transducer  is 
located  at  point  S.  Reflected  or  transmitted  waves  are  recorded  at  points  R  or 
T  Gross  sections  of  two  steps  used  are  shown  below. 


Figure  2.  Waveforms  of  Rayleigh  waves  transmitted  across  a  3mm  vertical  step. 
The  top  trace  is  the  input  signal.  The  small  pulse  arriving  1.0  psec  before  the 
Rayleigh  wave  is  an  S  wave  excited  at  the  step.  The  vertical  scales  on  the 
seismograms  have  been  expanded  by  the  scaling  factors  shown  at  the  start  oT 
each  trace.  A:  Transmission  in  the  up-step  direction.  B:  Transmission  in  l.h 
down-step  direction. 


Figure  3.  Energy  transmission  coefficients  for  the  Rayleigh  waves  in  Figure  2. 

The  horizontal  axis  is  normalized  frequency  (/=—),  where  h  is  the  step  height. 

A 

and  A  is  the  wavelength.  Curves  are  shown  for  three  angles  of  incidence.  A: 
Transmission  in  the  up-step  direction.  B:  Transmission  in  t.he  down-step  direc¬ 
tion. 


Figure  4  Waveforms  of  Rayleigh  waves  reflected  from  the  3mrn  vertical  step. 
Note  the  secondary  reflection  corresponding  to  a  reflection  from  the  upper 
edge  of  the  step  A:  Reflection  from  the  up-step.  B:  Reflection  from  the  down- 
step. 


Figure  f).  Energy  reflection  coefficients  for  the  reflected  Rayleigh  waves  in  Fig¬ 
ure  4.  Note  t  he  change  of  vertical  scale  from  Figure  3.  A  Reflection  from  the 
up-step  B:  Reflection  from  the  down-step. 


Figure  6  Seismograms  of  Rayleigh  waves  transmitted  across  (A)  and  reflected 
from  (B)  th  e  '.inrun,  60°  ramp.  Propagation  is  in  the  up-step  direction  The 
seismograms  in  (B)  were  high-pass  Altered  at  0.3  MHz  at  the  time  of  recording 
to  I'liliaiiec  the  signal-t o-noise  ratio. 


Figure  7.  Energy  transmission  and  reflection  coefficients  for  the  Rayleigh  waves 
m  Figure  0.  A:  tip-step  transmission.  B:  Up-step  reflection.  Note  the  difference 
iri  t  ho  umpli  l.i  ide  v.f,t|es. 


figure  6.  Energy  t runsnnssioii  data  of  Figure  3a  after  application  of  the  empiri¬ 
cal  relation  i » f  Equation  (3).  The  t  /  ansmissioti  curves  Tor  the  different  incidence 
ii  r ,  v  1 » ■  <  MOW  v  < )  i  r  t  f  ‘ !  d  •  * 


Figure  0.  The  I  of  t  he  3rnm  vertical  step  (A)  and  ramp  (B)  models  on  a 

dispersed  wave  tram.  I'ht  input,  signal  (top  tr.u  e)  is  a  linear  chirped  sinusoid, 
who->  treijuenry  increases  from  0.2  to  2.0  Mliz  in  60  /usee.  The  corresponding 
r  angr  of  t  h>-  r  atio  of  step  height /wavelength  (/  )  is  0.2  to  2. 
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which  surface  waves  cross  regions  of  strong  structural  contrast.  In  basiri-ami- 
rarige  provinces,  for  example,  where  crustal  Rayleigh  waves  traverse  numerous 
steps,  strong  effects  can  be  expected  due  to  compounding  of  the  effects 
suffered  by  passage  through  a  single  step. 

The  ultrasonic  model-synthetic  method  may  prove  useful  in  more  general 
problems  of  three-dimensional  surface  wave  scattering,  particularly  in  situa¬ 
tions  for  which  analytic  or  numerical  solutions  are  difficult  to  obtain. 


model-synthetic  data. 

The  agreement  between  real  and  modeled  waveforms  is  generally  good, 
particularly  at  longer  periods  (Figure  13).  Note  that  the  relative  amplitudes 
between  SHI,  and  NDI  are  similar  for  both  the  observed  and  model-synthetic 
seismograms.  No  attempt  has  been  made  here  to  model  the  small-scale  scatter¬ 
ing  in  Tibet  arid  the  Himalayas;  the  short  period  (<15-20  sec)  components  of  the 
waveforms  therefore  do  not  match  as  well. 


CONCLUSIONS 

Using  three-dimensional  ultrasonic  modeling  techniques,  transmission  and 
reflection  coefficients  as  a  function  of  wavelength  and  incidence  angle  are 
obtained  for  Rayleigh  waves  propagating  in  simple  step  models.  The  efficiency  of 
Rayleigh  wave  transmission  and  reflection  varies  strongly  with  both  incidence 
angle  arid  wavelength.  At  some  incidence  angles  and  frequencies  as  much  as 
90%  of  the  incident  Rayleigh  wave  energy  is  scattered  into  body  waves. 
Transmission  coefficients  obtained  at  normal  (90°)  incidence  agree  with  those 
obtained  with  two  dimensional  modeling  methods  and  two-dimensional  finite- 
difference  methods.  An  empirical  relation  between  transmission  coefficients  at 
different,  incidence  angles  is  obtained  for  the  case  of  a  vertical  step.  This  rela¬ 
tion  may  be  applied  to  two-dimensional  theoretical  (numerical)  calculations  of 
propagation  in  more  complicated  structures  to  obtain  estimates  of  Rayleigh 
wave  transmission  coefficients  at.  oblique  incidence  angles. 

The  effect  of  simple  structures  on  Rayleigh  wave  propagation  in  the  Earth 
is  demonstrated  by  convolving  the  impulse  response  of  the  models  with  real  or 
synthetic  dispersed  wave  trains.  Some  observed  features  of  Rayleigh  waves 
crossing  an  ocean-continent  margin  and  the  Tibetan  Plateau  are  similar  to 
those  predicted  from  model  results.  The  method  may  be  applied  to  studies  in 
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angle  governed  by  Snell’s  law  (e.g.  Evernden,  1953;  Capon,  1970).  Using  the 
coast  and  ocean  models  given  by  Drake  and  Dolt  (1980)  we  find  that,  in  the 
period  range  of  interest  in  this  example  (15-40  sec),  Rayleigh  waves  incident  on 
the  coast,  at  incidence  angles  greater  than  30°  are  refracted  laterally  by  no 
more  than  about  8°.  This  is  within  the  uncertainty  (fl0°)  in  the  estimate  of 
incidence  angle  and  so  is  not  an  important  consideration  in  this  example. 

THE  TiBh'TAN  PLATKAU 

The  Tibetan  Plateau  is  often  regarded  as  an  isostatieally-compensated 
block  of  uplifted  crustal  rock  with  an  average  elevation  of  5  km  (e.g.  Bird  and 
Tokso/,  1977).  According  to  the  ultrasonic  modeling  results,  we  expect  that 
Rayleigh  waves  crossing  the  margins  of  such  a  block  will  be  scattered  at  the 
edges  >f  the  block  and  the  waveforms  distorted.  To  investigate  the  effect  of  the 
Tibet.a  i  Plateau  on  Rayleigh  waves  propagating  across  it  we  compare  real  and 
model- synthetic  Rayleigh  wave  seismograms  for  the  propagation  paths  shown  in 
Figure  12. 

Ti  e  Tibetan  Plateau  is  modeled  here  as  a  cylindrical  lens  with  effective 
height  of  10km  above  average  terrain  (Figure  12,  inset).  The  lens  itself  lies  on  a 
half-spare.  The  lens  includes  the  effects  of  both  the  elevated  topography  of  the 
Tibet. a  i  Plateau  arid  the  deep  crustal  root  of  the  Himalayas.  Model-synthetic 
Rayleii.h  seismograms  are  generated  by  convolving  an  input  waveform  with  the 
impels  ■  r  esponse  of  the  step  model  appropriate  for  each  of  the  two  edges  of  the 
fens.  l'-e.iuse  the  Tibetan  Plateau  is  wide  compared  t.o  a  Rayleigh  wavelength, 
mt'Tf'  renre  arising  from  multiple  reflections  from  opposite  edges  of  the  Pla- 
''  a  u  i  m  gleet  ed  Since  there  are  no  WWSSN  stations  on  the  northern  front  of 
the  PI  trail,  we  use  the  Rayleigh  wave  recorded  at  KDI,  as  an  input  signal  t.o  the 
1‘  n  ■  ii  mlel  Hu  regional  at  ten  nut  ion  model  of  Bird  (1978)  is  applied  to  the 


THE  OCEAN-  CONTINENT  MARGIN 

The  incidence  angle  dependence  of  Rayleigh  wave  transmission  across  the 
ocean-continent,  margin  can  be  seen  in  recordings  of  teleseismic  oceanic  Ray¬ 
leigh  waves  made  at  coastal  stations.  For  example,  Figure  11  shows  Rayleigh 
waves  recorded  at  W'WSSN  stations  on  the  west  coast  of  the  United  States  from 
an  earthquake  in  the  Loyalty  Islands.  The  incidence  angle  is  taken  to  be  the 
angle  between  the  great-circle  path  and  the  continental  shelf.  For  paths  l  > 
BKS,  COR  and  LON,  the  angles  are  90°,  60°,  and  45°,  respectively.  Among  the 
continental  stations  in  this  example,  epicentral  distance  ranges  from  87.2" 
(BKS)  to  91.9"  (LON).  Great-circle  azimuths  vary  by  no  more  than  about  1"  and 
source  radiation  patterns  would  have  the  same  effect  at  these  stations. 
Differences  between  the  seismograms  recorded  at  the  continental  stations  are 
therefore  primarily  due  to  variations  in  the  incidence  angle.  At  normal 
incidence  (BKS)  most  of  the  Rayleigh  wave  energy  is  confined  to  the  first  dozen 
or  so  cycles,  after  which  time  the  amplitude  decays  rapidly  to  a  small  value.  As 
the  incidence  angle  decreases  to  near-grazing,  the  initial  ’packet'  of  energy 
broadens  and  the  rate  of  amplitude  decay  decreases.  This  trend  is  consistent 
with  that  observed  in  the  chirped  sine  synthetics  (Figure  9). 

Synthetic  Rayleigh  wave  seismograms  were  generated  for  a  pure  oceanic 
path  using  the  PFJM-0  model  (Dziewonski  et  a l  ,  1975)  and  convolved  with  the 
impulse  response  of  the  3mm  ramp  models,  scaled  to  a  height  of  60  km  (Figure 
lit)).  The  model-synt hetic  seismograms  exhibit  azimuthal  behavior  similar  to 
that,  of  the  real  seismograms  (Figure  1  la). 

In  this  example  we  have  assumed  the  Rayleigh  waves  propagate  on  the 
great  circle  path  between  source  and  receiver  In  fact,  the  velocity  contrast 
across  the  ocean-cord  .’merit,  margin  in  the  Kart.h  is  such  that  Rayleigh  waves  of  a 
given  period  traversing  the  margin  will  deviate  from  the  great-circle  path  by  an 
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and  Tx,  /-a-,  and  Vx  are  characteristic  times,  lengths,  and  velocities,  respectively 


in  either  ttie  model  ( m )  or  the  Earth  (e). 


In  general,  the  terrestrial  scale  length  appropriate  for  a  given  physical 
problem  is  considerably  greater  than  that  which  would  be  indicated  by  surface 
topography  alone.  In  discussions  of  terrestrial  scaling,  the  terrestrial  scale 
length  is  not  to  be  interpreted  literally  as  the  dimension  of  any  real  structure; 
it  is  simply  the  characteristic  dimension  of  a  highly  simplified  model.  The  topo¬ 
graphic  expression  of  the  physical  model  includes  the  combined  effects  of  ter¬ 
restrial  topography  and  deep  structure.  The  ultrasonic  step  models  described 
in  the  first  section  of  this  paper  represent  a  very  simplified  approximation  to 
complex  boundaries  of  vertical  heterogeneities  of  the  earth  structure.  It  is  still 
worthwhile,  however,  to  investigate  whether  models  produce  effects  similar  to 
those  observed  in  the  earth. 


"We  fake  as  the  input  signal  for  a  given  model  a  real  or  synthetic  terrestrial 
seismogram.  The  seismogram  is  resampled  at  a  rate  determined  by  the  model 
scaling  factor,  and  convolved  with  the  impulse  of  the  model.  The  result  of  the 
convolution  is  resampled  to  restore  the  seismogram  to  its  original  (terrestrial) 
time  scale. 

Because  the  usable  frequency  of  Rayleigh  waves  in  the  modeling  experi¬ 
ments  is  limited  to  the  band  0.2— 2.0  MHz,  the  range  of  normalized  frequencies 
available  for  scale  modeling  is  also  limited.  The  range  of  step  heights  and 
wavelengths  which  may  be  modeled  reliably  given  the  band-limited  nature  of 


the  ultrasonic  pulse  is  shown  in  Figure  10. 


impulse  response  of  a  given  model  for  waves  propagating  in  a  given  direction  is 
thus  the  result  of  deconvolving  the  transducer/ recording  system  impulse 
response  from  the  recorded  signal.  In  the  frequency  domain  we  may  write 
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wtiere  fm(0;cj)  is  the  complex  transfer  function  for  the  desired  model  at 
incidence  angle  0,  Sm(Q\cj)  is  the  complex  Fourier  spectrum  of  a  Rayleigh  wav 
propagating  across  the  model  at.  angle  0,  5>.  (sj)  is  the  half-space  Rayleigh  wave 
spectrum,  and  v  is  the  angular  frequency.  The  inverse  Fourier  transform  of 
this  quantity  yields  the  impulse  response  of  the  given  model. 


As  a  test,  of  the  effect  of  a  step  on  a  dispersed  wave  train,  a  chirped 
sinusoid  was  convolved  with  the  impulse  response  of  the  3mrn  vertical  step  and 
ramp  models  at.  several  incidence  angles  (Figure  9).  At  some  incidence  angles 
and  frequencies  the  envelope  of  the  wave  train  drops  to  near  zero,  suggestive  of 
the  "beating”  phenomenon  frequently  observed  in  seismograms  of  Rayleigh 
waves  propagating  aeross  ocean  basins.  This  beating  is  usually  attributed  to 
mu  itipa  thing  due  to  lateral  heterogeneities  along  the  propagation  path 
(Kvernden  1 9b'd,  1 9b4).  At  least  part  of  the  amplitude  modulation  of  dispersed 
Rayleigh  wave  trains  observed  in  the  Earth  may  be  due  to  propagation  across 
sharp  structural  features  such  as  those  modeled  here. 

In  order  to  investigate  the  effect  of  simple  steps  on  Rayleigh  wave  propaga¬ 
tion  m  the  Earth,  the  physical  dimensions  of  the  model  are  sealed  to 
correspond  with  those  of  some  realistic  Earth  structure.  The  scaling  factor 
relating  Ihr  physical  di rriensions  of  t  tie  Earth  to  those  of  the  model  is  given  hy 
th<  linear  relation  (White,  19f>f>): 


Tm.  -  k  Te 


where  the  scaling  factor  k  is  given  by 


(5) 


using  the  empirical  r  ction 


/min(Q)  = 


/  mm(90°) 
sm0 
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where  /mjn(0)  is  the  value  of  f  at  which  the  first  minimum  in  the  transmission 
curve  appears.  When  this  relation  is  applied  to  the  entire  transmission  curve 
the  curves  at  each  incidence  angle  can  be  made  to  coincide  (Figure  B).  In 
terms  of  the  effective  step  height  at  angle  0,  Equation  (£a)  implies 


h.aff  (0)  =  h  sin0  (3b) 

To  a  first  approximation  the  transmission  coefficient  for  an  arbitrary  incidence 

angle  can  thus  be  estimated  from  the  normal  incidence  transmission  coefficient 

using  this  empirical  relation.  This  approximation  is  poor,  however,  at  incidence 

angles  shallower  than  about  30". 


Because  of  its  strong  incidence  angle  dependence,  the  reflection 
coefficient  for  the  step  is  difficult  to  model  with  a  simple  empirical  relation. 


IMPLICATIONS  TOR  THE  EARTH 

One  of  the  most  interesting  features  of  the  ultrasonic  results  is  that  both 
transmitted  and  reflected  Rayleigh  wave  pulse  shapes  are  severely  altered  as  a 
result  of  interaction  with  the  various  steps.  The  waveform  distortion  is  due  to 
the  combined  effects  of  frequency-dependent  energy  transmission  and 
reflection  and  interference  of  waves  multiply  reflected  within  the  step.  More 
complicated  input,  signals,  such  as  dispersed  wave  trains,  can  be  expected  to 
undergo  similar  kinds  of  angle-dependent  waveform  distortion  after  propagat¬ 
ing  across  ,i  step. 

The  effect  of  a  given  model  on  a  surface  wave  signal  is  determined  by  con¬ 
volving  the  desired  input,  signal  with  the  impulse  response  of  the  model.  As 
noted  earlier,  the  Rayleigh  wave  propagating  across  the  aluminum  half-space 
represents  the  impulse  response  of  the  transducers  and  recording  system.  The 
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Note  that  in  the  ease  of  the  rarnp  model  the  frequencies  have  been  normal¬ 
ized  to  a  step  height  of  A  =3.0  mm;  in  fact  this  normalization  is  arbitrary.  While 
the  vertical  change  in  elevation  is  indeed  3.0  mm,  the  distance  from  the  lower 
edge  to  the  upper  edge  as  measured  along  the  rarnp  surface  is  actually  larger 

than  this  (——• — ,  where  a  is  the  ramp  angle).  This  distance  may  also  enter  into 
sma 

the  scaling  relations.  Direct  comparisons  between  step  and  ramp  transmission 
and  reflection  coeffcicnts  may  not,  therefore,  be  appropriate  when  scaled  t.c 
the  same  stop  height. 

AZIMUTHAL  DEPENDENCE  OF  SCALING 

In  the  models  discussed  above,  the  character  of  the  transmitted  and 
reflected  waveforms  changes  most  rapidly  as  the  incidence  angle  decreases 
below  about  30°.  It  is  in  the  range  0°  to  30°  that  the  behavior  of  Rayleigh  wave 
propagation  departs  most  strikingly  from  that  predicted  by  two-dimensional 
theoretical  and  approximate  methods.  At  incidence  angles  greater  than  SO3 
simple  empirical  relations  can  be  found  between  the  results  obtained  at  normal 
and  oblique  incidence. 

In  the  case  of  transmission  through  the  vertical  step  (Figure  3),  note  that 
the  energy  transmission  minimum  near  /  =0.5  at  normal  incidence  moves  to 
higher  frequency  as  the  incidence  angle  decreases.  The  equivalent  effect  is 
observed  for  normally  incident  waves  when  the  step  height  is  decreased.  This 
suggests  that  the  effective  step  height  is  dependent  on  incidence  angle.  That  is, 
as  the  incidence  angle  decreases  the  incident  Rayleigh  wave  ''sees''  a  shorter 
step.  The  functional  dependence  of  effective  step  height  on  incidence  angle  is 
estimated  by  finding  a  simple  functional  relation  between  incidence  angle  and  a 
particular  feature  in  the  transmission  curves.  The  normalized  frequency  at 
which  the  first,  minimum  occurs  in  the  transmission  curves  is  approximated  well 
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